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1 The 1-D Heat Equation 
1.1 Physical derivation 

Reference: Guenther & Lee §1.3-1.4, Myint-U & Debnath §2.1 and §2.5 
[Sept. 8, 2006] 

In a metal rod with non-uniform temperature, heat (thermal energy) is transferred 
from regions of higher temperature to regions of lower temperature. Three physical 
principles arc used here. 

1. Heat (or thermal) energy of a body with uniform properties: 

Heat energy = emu, 

where m is the body mass, u is the temperature, c is the specific heat, units [c] = 
j^2rp-2jj~i units are M mass, L length, T time, U temperature), c is the energy 

required to raise a unit mass of the substance 1 unit in temperature. 

2. Fourier's law of heat transfer: rate of heat transfer proportional to negative 
temperature gradient. 

Rate of heat transfer "^it 

= -Ko^ (1) 

area ox 

where Kq is the thermal conductivity, units [Kq] = MLT~^U~^. In other words, heat 
is transferred from areas of high temp to low temp. 

3. Conservation of energy. 

Consider a uniform rod of length I with non-uniform temperature lying on the 
X-axis from x — to x — I. By uniform rod, we mean the density p, specific heat c, 
thermal conductivity Kq, cross-sectional area A are ALL constant. Assume the sides 
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of the rod are insulated and only the ends may be exposed. Also assume there is no 
heat source within the rod. Consider an arbitrary thin slice of the rod of width Ax 
between x and x + Ax. The slice is so thin that the temperature throughout the slice 
is u {x, t). Thus, 

Heat energy of segment = c x pAAx x u = cpAAxu {x,t) . 

By conservation of energy, 

change of 

heat in from heat out from 
heat energy of = — 

left boundary right boundary 

segment in time At 



Prom Fourier's Law (1), 

cpAAxu (x, t + At)- cpAAxu (x, t) = AtA ( -Ko^] - AtA ( -Ko^] 

\ dxj^ V 9xJ^^^^ 

Rearranging yields (recall p, c. A, Kq are constant), 

U{x,t + At) -U{x,t) (if):c+Aa: ~ 

At ~ cp \ Ax 

Taking the limit At, Ax — > gives the Heat Equation, 

du d'^u , , 

where 

cp 

is called the thermal diffusivity, units [k] = I? jT . Since the shce was chosen arbi- 
trarily, the Heat Equation (2) applies throughout the rod. 



1.2 Initial condition and boundary conditions 

To make use of the Heat Equation, we need more information: 

1. Initial Condition (IC): in this case, the initial temperature distribution in the 
rod u (x, 0). 

2. Boundary Conditions (BC): in this case, the temperature of the rod is affected 
by what happens at the ends, x — 0, Z. What happens to the temperature at the 
end of the rod must be specified. In reahty, the BCs can be complicated. Here we 
consider three simple cases for the boundary at x = 0. 
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(I) Temperature prescribed at a boundary. For t > 0, 



PDE: 


Uf — KUxxi 


< X <l. 


IC: 


u{x,^) = fix), 


< X < I, 


BC: 


u{0,t) =u{L,t) = 0, 


t > 0. 



u{0,t) = ui (t) . 

(II) Insulated boundary. The heat flow can be prescribed at the boundaries, 

-Ko^m^<l>r(t) 

(III) Mixed condition: an equation involving ^(O,^), du/dx {0,t), etc. 
Example 1. Consider a rod of length / with insulated sides is given an initial 

temperature distribution of / (x) degree C, for < x < /. Find u [x, t) at subsequent 
times t > if end of rod are kept at 0° C. 

The Heat Eqn and corresponding IC and BCs are thus 

(4) 
(5) 
(6) 

Physical intuition: we expect — > as i — > oo. 
1.3 Non-dimensionalization 

Dimensional (or physical) terms in the PDE (2): k, I, x, t, u. Others could be 
introduced in IC and BCs. To make the solution more meaningful and simpler, we 
group as many physical constants together as possible. Let the characteristic length, 
time and temperature be L*, and f/*, respectively, with dimensions [L*] = L, 
[T^] = T, = U . Introduce dimensionless variables via 

X - t -\ u(x.t) fix) 

t-TfT: U{x,t)=^-^, f(x) = '-^. (7) 

The variables x, i, u are dimensionless (i.e. no units, [x] — 1). The sensible choice 
for the characteristic length is L^ — l, the length of the rod. While x is in the range 
< X < Z, X is in the range < x < 1. 

The choice of dimensionless variables is an ART. Sometimes the statement of the 
problem gives hints: e.g. the length I of the rod (1 is nicer to deal with than I, an 
unspecified quantity). Often you have to solve the problem first, look at the solution, 
and try to simplify the notation. 
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Prom the chain rule, 



du du dt 
di^ *~didt 


f/* dii 




du ^ du dx 
dx * dx dx 


[/* du 


dx 


U.d^u 




Li dx"^ 





Substituting these into the Heat Eqn (4) gives 

du T^K d^u 

Ut = KUxx =^ 



dt LI dx^ 

To make the PDE simpler, we choose T* — Li/ k — / k, so that 

du d'^u 

< X < 1, t > 0. 



di dx''' 

The characteristic (diffusive) time scale in the problem is = /^/ac. For different 
substances, this gives time scale over which diffusion takes place in the problem. The 
IC (5) and BC (6) must also be non-dimensionalized: 

IC: -u (£, 0) = / (£) , 0<x<l, 
BC: u (0, i) = u (1, i) =0, i> 0. 

1.4 Dimensionless problem 

Dropping hats, we have the dimensionless problem 



PDE 
IC 
BC 



Ut = u^x, < a; < 1, (8) 

(x, 0) = / (x) , 0<x<l, (9) 

M(0,t) = = 0, t > 0, (10) 



where x, t are dimensionless scalings of physical position and time. 

2 Separation of variables 

Ref: Guenther & Lee, §4.2 and 5.1, Myint-U & Debnath §6.4 
[Sept 12, 2006] 

We look for a solution to the dimensionless Heat Equation (8) - (10) of the form 

u{x,t)= X{x)T{t) (11) 
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Take the relevant partial derivatives: 



^X"{x)T{t), Ut^X{x)T'{t) 



where primes denote differentiation of a single- variable function. The PDE (8), Ut — 
Uxx, becomes 

r (t) _ X" (x) 

T{t)~ X {x) 

The left hand side (l.h.s.) depends only on t and the right hand side (r.h.s.) only 
depends on x. Hence if t varies and x is held fixed, the r.h.s. is constant, and hence 
T' /T must also be constant, which we set to —A by convention: 

r (t) X" {x) 



T (t) X (x) 
The BCs become, for t > 0, 



—A, A — constant. (12) 



u (0, t) ^ X (0) T (i) = 

u{l,t) = X(l)T(t) = 

Taking T (t) — would give u — ior all time and space (called the trivial solution) , 
from (11), which does not satisfy the IC unless f (x) — 0. If you are lucky and 
/ (x) = 0, then u = is the solution (this has to do with uniqueness of the solution, 
which we'll come back to). If / {x) is not zero for all < x < 1, then T (t) cannot be 
zero and hence the above equations are only satisfied if 

X{0)=X{1) = 0. (13) 

2.1 Solving for X (x) 

Ref: Guenther & Lee, §4.2 and 5.1 and 7.1, Myint-U & Debnath §7.1 - 7.3 
We obtain a boundary value problem for X (x), from (12) and (13), 

X" (x) + XX{x) = 0, 0<x<l, (14) 

X(0)=X(1) = 0. (15) 

This is an example of a Sturm-Liouville problem (from your ODEs class). 
There are 3 cases: A > 0, A < and A = 0. 
(i) A < 0. Let A = -k^ < 0. Then the solution to (14) is 

X = Ae^^^ + Se-*^^ 
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for integration constants A, B found from imposing the BCs (15), 

X{0)^A + B = 0, X (1) = Ae'' + Be''' = 0. 

The first gives A — —B, the second then gives A [e^^ — l) =0, and since \k\ > we 
have A — B — u — Q, which is the trivial solution. Thus we discard the case A < 0. 

(ii) A = 0. Then X{x)^Ax + B and the BCs imply = X (0) = S, = X (1) = 
A, so that A — B — u — Q. We discard this case also. 

(iii) A > 0. In this case, (14) is the simple harmonic equation whose solution is 

X [x) — A cos + Bsm (VXx) . (16) 

The BCs imply = X (0) = A, and 5 sin = 0. Wc don't want 5 = 0, since that 
would give the trivial solution m = 0, so we must have 

sinVA = 0. (17) 

Thus = nvr, for any nonzero integer n {n — 1,2,3,...). We use subscripts 
to label the particular n-value. The values of A are called the eigenvalues of the 
Sturm- Liouville problem (14), 

A„ = n^TT^, n = 1,2,3,... 

and the corresponding solutions of (14) are called the eigenfunctions of the Sturm- Liouville 
problem (14), 

Xn{x) — bn sin (nTTx) , n— 1,2,3, ... (18) 
We have assumed that n > 0, since n <0 gives the same solution as n > 0. 

2.2 Solving for T (t) 

When solving for X [x), we found that non-trivial solutions arose for A = n^7r^ for all 
nonzero integers n. The equation for T {t) is thus, from (12), 

r (t) = -n\^T {t) 

and, for n, the solution is 

Tn = c„e-"''^'*, n = l,2,3,... (19) 
where the c^'s are constants of integration. 
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2.3 Full solution u{x,t) 

Ref: Myint-U & Debnath §6.4, Ch 5 

Putting things together, we have, from (11), (18) and (19), 

Wn(x,i) = S„sin(n7rx)e-"''''*, n = l,2,3,... (20) 

where Bn = Cnhn- Each function Un (a:, t) is a solution to the PDE (8) and the BCs 
(10). But, in general, they will not individually satisfy the IC (9), 

Un (x, 0) = Bn sin {m:x) — f (x) . 

We now apply the principle of superposition: if Ui and U2 are two solutions to the 
PDE (8) and BC (10), then CiUi + C2U2 is also a solution, for any constants Ci, C2. 
This relies on the linearity of the PDE and BCs. We will, of course, soon make this 
more precise.... 

Since each Un {x, 0) is a solution of the PDE, then the principle of superposition 
says any finite sum is also a solution. To solve the IC, we will probably need all the 
solutions Un, and form the infinite sum (convergence properties to be checked), 

oc 

u{x,t) = ^ u„(x,t) . (21) 

n=l 

u {x,t) satisfies the BCs (10) since each Un {x,t) does. Assuming term-by-term dif- 
ferentiation holds (to be checked) for the infinite sum, then u (x, t) also satisfies the 
PDE (8). To satisfy the IC, we need to find S„'s such that 

00 00 
f (x) ^ u{x,0) ^J2'^-ix,0) 

sin (mrx) . (22) 

n=l n=l 

This is the Fourier Sine Series of / (x) . 

To solve for the -B„'s, we use the orthogonality property for the eigenfunctions 
sin (mrx) , 

\ m n 1 

sm [mnx) sm [nnx) ax = < = -Omn (23) 

1/2 m — n 2 



where 5mn is the kronecker delta, 

8n 



m ^ n 

1 m = n 



The orthogonahty relation (23) is derived by substituting 

2 sin {rmix) sin {rmx) = cos ((m — n) nx) — cos ((m + n) nx) 
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into the integral on the left hand side of (23) and noting 

cos {rmrx) dx — Smo- 



I 

Jo 
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The orthogonality of the functions sin {nnx) is analogous to that of the unit vectors 
X and y in 2-space; the integral from to 1 in (23) above is analogous to the dot 
product in 2-space. 

To solve for the S„'s, we multiply both sides of (22) by sin {rmrx) and integrate 
from to 1: 

„1 oo „i 

/ sin (mnx) f (x) dx = -B„ / sin (nnx) sin {mnx) dx 
Jo „=i Jo 

Substituting (23) into the right hand side yields 

sin (mTTx) f (x) dx = Bn-5r, 



/ 

Jo 



1 2 -'nm 

n=l 



By definition of Snm, the only term that is non-zero in the infinite sum is the one 
where n — m, thus 



1 1 

sin {mnx) f {x) dx = -Bm 
2 



Rearranging yields 

B. 



2 sin (mTTx) f {x) dx. (24) 
Jo 

The full solution is, from (20) and (21), 

oo 

u (x, t)^J2^n sin (nnx) e""''''*, (25) 

n=l 

where Bn are given by (24). 

To derive the solution (25) of the Heat Equation (8) and corresponding BCs 
(10) and IC (9), we used properties of linear operators and infinite series that need 
justification. 



3 Example : Cooling of a rod from a constant ini- 
tial temperature 

Suppose the initial temperature distribution / (x) in the rod is constant, i.e. / (x) = 
Uq. The solution for the temperature in the rod is (25), 

oo 

u {x, t) = ^2 sin (mrx) e~" ^ *, 

n=l 
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where, from (24), the Fourier coefficients are given by 

Bn — 2 sin (nnx) f (x) dx — 2uo / sin (nnx) dx. 
Jo Jo 

Calculating the integrals gives 

D o ■ ( w o cos(n7r)-l 2uq ^ J 

Bn = 2uo / sm [mrx) dx = —2uo ^ — = ~ ^> = i 

Jo "-TT nn \^ 

In other words, 

B2n = 0, B2n-1 



n even 
^ nodd 



(2n - 1) TT 
and the solution becomes 

, , 4mo v-^ sin ((2n — 1) TTx) / .,29\ /^^x 

^ ^) = V ^ (2n - 1 ^""P ^^"^ ~ ^) • ^^^^ 

n=l ^ 

3.1 Approximate form of solution 

The number of terms of the series (26) needed to get a good approximation for u {x, t) 
depends on how close t is to 0. The series is 

\ — TT^t I sin (Sttx) _9^2j 



u(x,t)^ — ( sin (nx) e"'^ * + ^ e'^^ * + 

71 \ 3 

The ratio of the first and second terms is 

I second term! e"^'^^* |sin37rx| 



first term 3 sinTrx 



—Hit f I III 

< e using |sinnx| < n |sinx| 

< for t > 



< 0.00034 



It can also be shown that the first term dominates the sum of the rest of the terms, 
and hence 

u (x, t)^— sin (ttx) e-^'\ for t>^. (27) 

What does t — I/tt^ correspond to in physical time? In physical time t' — IH/k 
(recall our scaling - here we use t as the dimensionless time and t' as dimensional 
time), t — l/n^ corresponds to: 

t' ~ 15 minutes, for aim rod of copper (k ^ 1.1 cm^ sec~^) 
t' ~ 169 minutes, for a 1 m rod of steel {k, ^ 0.1 cm^ sec~^) 
t' ~ 47 hours, for a 1 m rod of glass (k ~ 0.006 cm^ sec~^) 
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At i = I/tt^, the temperature at the center of the rod {x — 1/2) is, from (27), 

4 

u {x, t) a; — Uq — 0.47^0. 
ne 

Thus, after a (scaled) time t — I/tt^, the temperature has decreased by a factor of 
0.47 from the initial temperature Uq. 

3.2 Geometrical visualization of the solution 

[Sept 14, 2006] 

To analyze qualitative features of the solution, we draw various types of curves in 

2D: 

1. Spatial temperature profile, given by 

u — u{x, to) 

where to is a fixed value of x. These profiles are curves in the ux-plane. 

2. Temperature profiles in time, 

u — u (xo, t) 

where xo is a fixed value of x. These profiles are curves in the lii-plane. 

3. Curves of constant temperature in the xi- plane (level curves), 

u (a;, t) = C 

where C is a constant. 

Note that the solution u = u{x, t) is a 2D surface in the 3D wxt-space. The above 
families of curves are the different cross-sections of this solution surface. Drawing the 
2D cross-sections is much simpler than drawing the 3D solution surface. 

Sketch typical curves: when sketching the curves in 1-3 above, we draw a few 
typical curves and any special cases. While math packages such as Matlab can be 
used to compute the curves from, say, 20 terms in the full power series solution (26), 
the emphasis in this course is to use simple considerations to get a rough idea of what 
the solution looks like. For example, one can use the first term approximation (27), 
simple physical considerations on heat transfer, and the fact that the solution u {x, t) 
is continuous in x and t, so that if ti is close to ti, u (x, ti) is close to u (x, ^2)- 

3.2.1 Spatial temperature profiles 

For fixed t — to, the first term approximate solution (27) is 

u {x, t) ^ ^6-'^'*° sin {ttx) , t > I/tt^. (28) 
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Figure 1: Spatial temperature profiles ^(a;, to)- 
This suggests the center of the rod, x = 1/2, is a hne of symmetry for u (x, t), i.e. 




and, for each fixed time, the location of the maximum (minimum) temperature if 
uq > {uq < 0). We can prove the symmetry property by noting that the original 
PDE/BC/IC problem is invariant under the transformation x — > 1 — x. Note also 
that, from (26), 

oo 

Ux {x, t) = Auq cos ((2n — 1) ttx) exp (— (2n — 1)^ vr^t) 

n=l 

Thus (1/2, t) = and Uxxi^f^jt) < 0, so the 2nd derivative test implies that 
a; = 1/2 is a local max. 

In Figure 1, we have plotted two typical profiles, one at early times t = to ~ and 
the other at late times t = to ^ 0, and two special profiles, the initial temperature 
at t — {u = Uq) and the temperature as t —>■ oo {u = 0). The profile for t = to ^ 
is found from the first term approximation (28). The line of symmetry x — 1/2 is 
plotted as a dashed line for reference in Figure 1. 
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Figure 2: Temperature profiles in time u {xQ,t). 

3.2.2 Temperature profiles in time 

Setting X — xo in the approximate solution (27), 

u{xo,t)Ki sin (ttxo)^ e"'''*, t>l/n^. 

Two typical profiles are sketched in Figure 2, one near the center of the rod {xq ~ 1/2) 
and one near the edges {xq pa or 1). To draw these we noted that the center of 
the rod cools more slowly than points near the ends. One special profile is plotted, 
namely the temperature at the rod ends {x = 0, 1). 

3.2.3 Curves of constant temperature (level curves of u{x,t)) 

In Figure 3, we have drawn three typical level curves and two special ones, u — 
(rod ends) and u — Uq (the initial condition). For a fixed x — Xq, the temperature 
in the rod decreases as t increases (motivated by the first term approximation (28)), 
as indicated by the points of intersection on the dashed line. The center of the rod 
(x = 1/2) is a line of symmetry, and at any time, the maximum temperature is at 
the center. Note that at t = 0, the temperature is discontinuous at ,t = 0, 1. 

To draw the level curves, it is easiest to already have drawn the spatial temperature 
profiles. Draw a few horizontal broken lines across your u vs. x plot. Suppose you 
draw a horizontal line u = Ui. Suppose this line u = Ui crosses one of your profiles 
u{x,tQ) at position x = Xi. Then (a;i,to) is a point on the level curve u{x,t) = Ui. 
Now plot this point in your level curve plot. By observing where the line u = Ui 
crosses your various spatial profiles, you fill in the level curve u {x,t) = ui. Repeat 
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u = 



u = C^O 



u = 




u = C ^ ^Uq 



u = C ^ Uq 



X 



Figure 3: Curves of constant temperature u{x,t) = c, i.e., the level curves of u{x,t). 

this process for a few values of ui to obtain a few representative level curves. Plot 
also the special cases: u {x,t) = uq/w, u = 0, etc. 

When drawing visualization curves, the following result is also helpful. 

3.2.4 Maximum Principle for the basic Heat Problem 

Ref: Guenther & Lee §5.2, Myint-U & Debnath §8.2 

This result is useful when plotting solutions: the extrema of the solution of the 
heat equation occurs on the space-time "boundary", i.e. the maximum of the initial 
condition and of the time-varying boundary conditions. More precisely, given the 
heat equation with some initial condition / (x) and BCs u{0,t), u{l,t), then on a 
given time interval [0, T], the solution u (x, t) is bounded by 



Ur 



< u (x, t) < Un 



where 



max < max f (x) , max u (0, t) , max u (1, t) 

' 0<x<l ■> ^ " o<t<T ^ ' ^ ' o<t<T ^ ^ 



Urain = miu < miu fix), min u(0,t), min u(l,t) 

' 0<x<l ^ ^ ' o<i<T ^ ' ^ ' 0<t<T ^ ' 
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Thus, in the example above, Umin — and x^max = Uq, hence for all x e [0, 1] and 
t e [0,r], 0<u{x,t)< uo. 

4 Equilibrium temperature profile (steady-state) 

Intuition tells us that if the ends of the rod are held at 0" C and there are no heat 
sources or sinks in the rod, the temperature in the rod will eventually reach 0. The 
solution above confirms this. However, we do not have to solve the full problem to 
determine the asymptotic or long-time behavior of the solution. 

Instead, the equilibrium or steady-state solution u — ue{x) must be independent 
of time, and will thus satisfy the PDE and BCs with ut — 0, 

u'l^ = 0, < x < 1; ue (0) = Ue (1) = 0. 

The solution is ue = Cix + C2 and imposing the BCs implies ue {x) = 0. In other 
words, regardless of the initial temperature distribution u {x, 0) = / (x) in the rod, 
the temperature eventually goes to zero. 

4.1 Rate of decay of u{x,t) 

How fast does u approach ue — 0? Prom our estimate (82) above, 

\un{x,t)\ < Be-"'^''^\ n = 1,2,3,... 

2 2 2 / 2 \ "' 

where S is a constant. Noting that e~'^ * < e~"^ * — i * ) , we have 

oo 



n=l 



< {x,t)\ <bY,t-^ r = e--'* <l{t> 0). 



n=l n=l 



The last step is the geometric series result Yl^=\f^ — iz^) for |r| < 1. Thus 



y^^n {X,t) 



n=l 



< . t > 0. (29) 



1 - e-'^^ 



Therefore 



• u{x,t) approaches the steady-state ue {x) = exponentially fast (i.e. the rod 
cools quickly) 

• the first term in the series, sin (ttx) e~'^^* (term with smallest eigenvalue X — tt) 
determines the rate of decay oi u{x, t) 

• the BnS may also affect the rate of approach to the steady-state, for other 
problems 
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4.2 Error of first-term approximation 



Using the method in the previous section, we can compute the error between the 
first-term and the full solution (26), 



\u {x, t) — Ui {x, t)\ 



y^Mn ix,t) -Ui{x,t) 



n=l 

oo 



n=2 



<5^K(a:,t)|<i?5^r"=- 



n=2 



n=2 



where r = e *. Hence the solution u {x, t) approaches the first term ui {x, t) expo- 
nentially fast. 



\u(x,t) — Ui {x,t)\ < 



1 ^ > 0- 



(30) 



With a little more work we can get a much tighter (i.e. better) upper bound. We 
consider the sum 

oo 



n=N 



Substituting for Un {x, t) gives 

oo oo oo 

J2 K (x,t)\ <BJ2 e-"''^'* = Be-""'^'' e'^"'"^')"'* 

n=N n=N n=N 

Here's the trick: for n> N, 

n^-N^ = (n + N) {n - N) > 2N {n - N) > 0. 
Since is a decreasing function, then 

g-(n2-Ar2),r2t ^ ^-2N{n-N)n'^t 

for t > 0. Using the inequahty (32) in Eq. (31) gives 

oo oo 



n=N 



n=N 



N 



-2NttH 



n=N 

Simplifying the expression on the right hand side yields 

2^ \Un{x,t)\ < : 
n=N 
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-2NttH 



(31) 



(32) 



(33) 



For N — 1, (33) gives the temperature decay 



\u{x,t) \ = 



Un {X, t) 



n=N 



(34) 



n=N 



which is shghtly tighter than (29). For N — 2, (33) gives the error between u{x,t) 
and the first term approximation, 



\u {x, t) — Ui {x, t)\ 



n=2 



n=2 



which is a much better upper bound than (30). 

Note that for an initial condition u (x, 0) = / [x) that is is symmetric with respect 
to X = 1/2 (e.g. / {x) = uq), then U2n — for all n, and hence error between u {x, t) 
and the first term ui {x, t) is even smaller, 



\u {x, t) — Ui {x, t)\ = 
Applying result (33) with N = 3 gives 



{X,t) 



n=3 



< ^\Un {X,t)\ 



n=3 



\u {x,t) — Ui {x,t)\ < 



Be 



1 - e-e^^'* ■ 



Thus, in this case, the error between the solution u {x, t) and the first term ui {x, t) 
decays as e"^'^^* - very quickly! 



5 Review of Fourier Series 

Ref: Guenther & Lee, §3.1, Myint-U & Debnath §5.1-5.3, 5.5-5.6 
[Sept 19, 2006] 

Motivation: Recall that the initial temperature distribution satisfies 

oo 

f {x) — u {x, 0) = Bn sin {mrx) . 

n=l 

In the example above with a constant initial temperature distribution, / {x) — uq, we 
have 

Auo^sm{{2n-l)7ix) 
= 2^1 ■ (^^^ 

n=l 

Note that at a; = and x = 1, the r.h.s. does NOT converge to Uq ^ 0, but rather to 
(the BCs). Note that the Fourier Sine Series of / (x) is odd and 2-periodic in space 
and converges to the odd periodic extension of / (x) = uq. 
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Odd periodic extension The odd periodic extension of a function / (x) defined 
for X e [0, 1], is the Fourier Sine Series of / (x) evaluated at any x e R, 

oo 

/ ( x) = ^ Bn sin (nnx) . 

n=l 

Note that since sin (nTrx) = — sin(— nTrx) and sm{mrx) — sin (n7r (x + 2)), then 
f (x) — —f{—x) and f{x + 2) — f (x). Thus f (x) is odd, 2-periodic and f (x) 
equals / (x) on the open interval (0, 1). What conditions are necessary for / (x) to 
equal / (x) on the closed interval [0, 1]? This is covered next. 

Aside: canceUing uq from both sides of (35) gives a really comphcated way of 
writing 1, 

4 sin {{2n — 1) irx) 
~ TT ^ 2n- 1 ■ 

n=l 

5.1 Fourier Sine Series 

Given an integrable function / (x) on [0, 1] , the Fourier Sine Series of / (x) is 

oo 

Bn sin (nnx) (36) 

n=l 

where ^ 

Bji — 2 f (x) sin (nnx) dx. (37) 
Jo 



The associated orthogonality properties are 



sin {rrnix) sin {mix) dx 



1/2, m^ny^O, 
0, m ^ n. 



5.2 Fourier Cosine Series 

Given an integrable function / (x) on [0, 1], the Fourier Cosine Series of / (x) is 



^0 + ^ COS {mrx) (38) 



n=l 



where 



^0 — fix)dx, (average of / (x)) (39) 
Jo 

An ^ 2 f {x)cos{mrx)dx, n>l. (40) 
Jo 
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The associated orthogonality properties are 
cos {rmrx) cos (mrx) dx - 



I 

Jo 



1/2, m = n^O, 
0, m ^ 11. 



e.g. The cosine series of / {x) — uq for x e [0, 1] is just uq. In other words, Aq — uq 
and An — Oioicn> 1. 

5.3 The (full) Fourier Series 

The Full Fourier Series of an integrable function / (x), now defined on [—1,1], is 

oo 

/ (x) = Oo + (a„ cos {nnx) + bn sin (nnx)) (41) 



n=l 



where 



Oo = ^ y f{x)dx 

f {x) cos {nnx) dx 

bn — J f (x) sin {nnx) dx 

The associated orthogonality properties of sin and cos on [—1, 1] are, for any m,n 
1,2,3,... 

J sin {mnx) cos {nnx) dx — 0, all m, n, 



J sin {mnx) sin {nnx) dx ■ 
J cos {mnx) cos (riTrx) 



1, m — n ^ 0, 

0, m ^ n, 

1, m = n 7^ 0, 
0, m ^ n. 



5.4 Piecewise Smooth 

Ref: Gucnthcr & Lcc p. 50, 

Provided a function / {x) is integrable, its Fourier coefficients can be calculated. 
It does not follow, however, that the corresponding Fourier Series (Sine, Cosine or 
Full) converges or has the sum / {x). In order to ensure this, / {x) must satisfy some 
stronger conditions. 
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Definition Piecewise Smooth function: A function / (x) defined on a closed 
interval [a, b] is said to be piecewise smooth on [a, b] if there is a partition of [a,b], 



a = Xq < Xi < X2 < ■ ■ ■ < X, 



n 



b 



such that / has a continuous derivative (i.e. C^) on each closed subinterval [x^, a^m+i]- 
E.g. any function that is on an interval [a, b] is, of course, piecewise smooth on 



is piecewise smooth on [0, 1] , but is not continuous on [0, 1] . 

E.g. the function / [x) — \x\ is both continuous and piecewise smooth on [—1, 1], 
despite /' {x) not being defined at x = 0. This is because we partition [—1, 1] into 
two subintervals [—1, 0] and [0, 1]. When worrying about /' {x) near x = 0, note that 
on [0, 1] we only care about the left limit /' (0~) and for [—1, 0], we only care about 
the right limit /'(0+). 

1 /2 

E.g. the function f [x) = \x\ is continuous on [—1, 1] but not piecewise smooth 
on [—1, 1], since /' (0~) and /' (0+) do not exist. 

5.5 Convergence of Fourier Series 

Ref: Guenther & Lee p. 49 and (optional) §3.3, Myint-U & Debnath §5.10 

Theorem [Convergence of the Fourier Sine and Cosine Series]: If / (x) is piecewise 
smooth on the closed interval [0, 1] and continuous on the open interval (0, 1), then 
the Fourier Sine and Cosine Series converge for all x e [0, 1] and have the sum / [x] 
for all X e (0,1). 

Note: Suppose / {x) is piecewise smooth on [0, 1] and is continuous on (0, 1) except 
at a jump discontinuity &X, x — a. Then the Fourier Sine and Cosine Series converge 
to / {x) on (0, 1) and converge to the average of the left and right hmits x — a, 
i.e. (/ (a—) + / (a+)) /2. At the endpoints x = 0, 1, the Sine series converges to zero, 
since sin {m:x) = at x = 0, 1 for all n. 

Theorem [Convergence of the Full Fourier Series] : If / {x) is piecewise smooth on 
the closed interval [—1, 1] and continuous on the open interval (—1, 1), then the Full 
Fourier Series converges for all x e [—1, 1] and has the sum / [x) for all x e (—1, 1). 

5.6 Comments 

Ref: see problems Guenther & Lee p. 53 



[a A- 



E.g. the function 




2x, < X < 1/2 
1/2, l/2<x<l 



19 



Given a function / (x) that is piecewise smooth on [0, 1] and continuous on (0, 1), 
the Fourier Sine and Cosine Series of / (x) converge on [0, 1] and equal / (x) on the 
open interval (0, 1) (i.e. perhaps excluding the endpoints). Thus, for any x e (0, 1), 

oo oc 

Aq + An COS {mrx) — f {x) — Bn sin {mrx) . 

n=l n=l 

In other words, the Fourier Cosine Series (left hand side) and the Fourier Sine Series 
(right hand side) are two different representations for the same function / (x), on the 
open interval (0, 1). The values at the endpoints a; = 0, 1 may not be the same. The 
choice of Sine or Cosine series is determined from the type of eigenfunctions that give 
solutions to the Heat Equation and BCs. 



6 Well-Posed Problems 

Ref: Guenther & Lee §1.8 and §5.2 (in particular p. 160 and preceding pages, and p. 
150) 

Ref: Myint-U & Debnath §1.2, §6.5 

We call a mathematical model or equation or problem well-posed if it satisfies the 
following 3 conditions: 

1. [Existence of a solution]: The mathematical model has at least 1 solution. 
Physical interpretation: the system exists over at least some finite time interval. 

2. [Uniqueness of solution]: The mathematical model has at most 1 solution. Phys- 
ical interpretation: identical initial states of the system lead to the same out- 
come. 

3. [Continuous dependence on parameters]: The solution of the mathematical 
model depends continuously on initial conditions and parameters. Physical 
interpretation: small changes in initial states (or parameters) of the system 
produce small changes in the outcome. 

If an IVP (initial value problem) or BIVP (boundary initial value problem - e.g. 
Heat Problem) satisfies 1, 2, 3 then it is well-posed. 

Example: for the basic Heat Problem, we showed 1 by construction a solution 
using the method of separation of variables. Continuous dependence is more difficult 
to show (need to know about norms), but it is true, and we will use this fact when 
sketching solutions. Also, when drawing level curves u {x, t) = const, small changes 
in parameters {x, t) leads to a small change in u. We now prove the 2nd part of 
well-posedness, uniqueness of solution, for the basic heat problem. 
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6.1 Uniqueness of solution to the Heat Problem 

Definition We define tlie two space-time sets 



V = {{x,t) : < X < 1, i > 0} 

V = {{x,t) : < X < 1, t > 0} 

and tlie space of functions 

(D) = [u{x,t) : Uxx continuous in V and u continuous in P j . 

In otfier words, tfie space of functions that are twice continuously differentiable on 
[0, 1] for i > and continuous on [0, 1] for t >0. 

Theorem The basic Heat Problem, i.e. the Heat Equation (8) with BC (10) 
and IC (9), 



PDE 
IC 
BC 



Ut = Uxx, < X < 1, 

M (x, 0) = / (x) , < X < 1, 

w(0,i) = 1^(1,^) = 0, t > 0, 



has at most one solution in the space of functions {V) . 

Proof: Consider two solutions -Ui, U2 E (X') to the Heat Problem. Let v = 
ui — U2- We aim to show that t; = on [0, 1], which would prove that ui = U2 and 
then solution to the Heat Equation (8) with BC (10) and IC (9) is unique. Since each 
of Ui, U2 satisfies (8), (9), and (10), the function v satisfies 

vt = (ui - U2)t (42) 

= Uit - U2t 
= {Ui-U2)^^ 

= Vxx, < a; < 1, 

and similarly, 

IC: V {x, 0) = Ml {x, 0) - U2 {x, 0) = / (x) - / (x) = 0, < x < 1, (43) 
BC: v{0,t) ^ui{0,t) -U2{0,t) ^0, v{l,t)^0, t > 0. (44) 



Define the function 



V{t)^ I (x, t)dx>^, t> 0. 
Jo 
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Showing that v {x, t) — reduces to showing that V (t) — 0, since v {x, t) is continuous 
on [0, 1] for ell t > and if there was a point x such that v {x, t) ^ 0, then V (t) 
would be strictly greater than 0. To show V (t) — 0, we differentiate V (t) in time 
and substitute for vt from the PDE (42), 



dV 

~dt 



! / vvtdx = 2 / 
Jo Jo 



VVxxdx 



Integrating by parts (note: vVxx = {'Wx)x ~ '^D gives 
dV 



^0 



Using the BCs (44) gives 







dy_ 

~dt 



-2 / vldx < 0. 



X 



V{0)= / v^x,0)dx = 0. 



The IC (43) implies that 



Thus, V{t) > 0, dV/dt < 0, and V (0) = 0, i.e. V (t) is a non-negative, non- 
increasing function of time whose initial value is zero. Thus, for all time, V (t) — 
and V {x, t) — for all x e [0, 1], implying that Ui — U2- This proves that the solution 
to the Heat Equation (8), its IC (9), and BCs (10) is unique, i.e. there is at most one 
solution. Uniqueness proofs for other types of BCs follows in a similar manner. 



7 Variations on the basic Heat Problem 

[Sept 21, 2006] 

We now consider variations to the basic Heat Problem, including different types 
of boundary conditions and the presence of sources and sinks. 



7.1 Boundary conditions 

7.1.1 Type I BCs (Dirichlet conditions) 

Ref: Guenther & Lee p. 149 

Type I, or Dirichlet, BCs specify the temperature u {x, t) at the end points of the 
rod, for t > 0, 

u{0,t) = gi{t), 
u{l,t) = g^it). 
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Type I Homogeneous BCs are 

u{0,t) = 0, 
u{l,t) = 0. 

The physical significance of these BCs for the rod is that the ends are kept at 0° C. 
The solution to the Heat Equation with Type I BCs was considered in class. After 
separation of variables u {x, t) — X {x) T {t), the associated Sturm- Liouville Boundary 
Value Problem for X {x) is 

X" + AX = 0; X(0)^X (1) = 0. 

The eigenfunctions are Xn (x) = Bn sin {httx). 

7.1.2 Type II BCs (Newmann conditions) 

Ref: Guenther & Lee p. 152 problem 1 

Type II, or Newmann, BCs specify the rate of change of temperature du/dx (or 
heat flux) at the ends of the rod, for t > 0, 

du ,^ . , , 

^(0,^) = 91 it), 

9u ,^ . 

Type II Homogeneous BCs are 

g(0,.) ^ 0, 

^ 0. 

The physical significance of these BCs for the rod is that the ends are insulated. 
These lead to another relatively simple solution involving a cosine series (see problem 
6 on PS 1). After separation of variables u {x, t) = X (x) T {t), the associated Sturm- 
Liouville Boundary Value Problem for X (x) is 

X" + XX = 0; X' (0) = X' (1) = 0. 

The eigenfunctions are Xq (x) — Aq — const and Xn (x) — An cos (mrx) . 

7.1.3 Type III BCs (Mixed) 

The general Type HI BCs are a mixture of Type I and II, for t > 0, 

OL\ ^ (0, t) + 0L2U (0, t) = gx {t) , 
aa^^ (l,t) + (l,t) = 92 (t). 
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After separation of variables u{x,t) — X{x)T{t), the associated Sturm-Liouville 
Boundary Value Problem for X (x) is 

X" + XX = 0, 

aiX' (0) + a2X (0) = 
asX' (0) + a4X (0) = 

The associated eigenfunctions depend on the values of the constants q;i,2,3,4- 
Example 1. ai = 0^4 = 1, a2 = 0:3 = 0. Then 

X" + AX = 0; X' (0) = X (1) = 

and the eigenfunctions are Xn = An cos (^^y^Tra;) . Note: the constant Xq = Aq is not 
an cigcnfunction here. 

Example 2. ai = 04 = 0, a2 = "3 = 1- Then 

X" + XX^O; X (0) = X' (1) = 

and the eigenfunctions are X„ = 5„sin (^^^ttx). Note: the constant Xq = Aq is not 
an cigcnfunction here. 

Note: Starting with the BCs in Example 1 and rotating the rod about x = 1/2, 
you'd get the BCs in Example 2. It is not surprising then that under the change of 
variables x = 1—x, Example 1 becomes Example 2, and vice versa. The eigenfunctions 
also possess this symmetry, since 

. /2n-l ,^ A , /2n-l \ 

sml^ — 7r(l -x) 1 = (-1) cosl^ — ttx j . 

Since we can absorb the (—1)" into the constant the eigenfunctions of Example 
1 become those of Example 2 under the transformation x — 1 — x, and vice versa. 

7.2 Solving the Heat Problem with Inhomogeneous (time- 
independent) BCs 

Ref: Gucnthcr & Lee p. 149 

Consider the Heat Problem with inhomogeneous Type I BCs, 

Ut — Uxx-i < X < 1 

ii(0,t) = 0, u{l,t) = ui = ccnst, i > 0, (45) 
M(a;,0) = 0, 0<a;<l. 
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Directly applying separation of variables u {x, t) — X (x) T (t) is not useful, because 
we'd obtain X (1) T (t) — ui for t > 0. The strategy is to rewrite the solution u {x, t) 
in terms of a new variable v{x,t) such that the new problem for v has homogeneous 
BCs! 

Step 1. Find the steady-state, or equilibrium solution ue{x), since this by defini- 
tion must satisfy the PDE and the BCs, 

u% ^ 0, 0<x <1 
ue (0) = 0, ue (1) = Ui — const. 

Solving for ue gives ue {x) = Uix. 

Step 2. Transform variables by introducing a new variable v {x,t), 

V {x, t) — u (x, t) — Ue {x) — u (x, t) — uix. (46) 

Substituting this into the Heat Problem (45) gives a new Heat Problem, 

Vt = Vxx, < X < 1 

v{0,t) = 0, v{l,t) = 0, t>0, (47) 
v{x,0) — —uix, < a; < 1. 

Notice that the BCs are now homogeneous, and the IC is now inhomogeneous. Notice 
also that we know how to solve this - since it's the basic Heat Problem! Based on our 
work, we know that the solution to (47) is 

oo „1 

V {x, t) = y2Bn sin {mix) e""^"^*, 3^ = 2 f (x) sin (mrx) dx (48) 

n=l -^0 

where / (x) = —Uix. Substituting for / (x) and integrating by parts, we find 



-2ui [ 
Jo 



Bn — —2ui / xsva.{m:x)dx 



X cos m^x 



nn 

2u, {-If 



x=0 JO 



COS (nnx) dx 



nn 

Step 3. Transform back to u{x,t), from (46), (48) and (49), 



(49) 



2u °° ( IV" 

u (x, t) = Ue (x) + V (x, t) = uix H — sin (mrx) e"'*^''^^ (50) 

TT ' n 

n=l 

The term ue (x) is the steady state and the term v {x, t) is called the transient, since 
it exists initially to satisfy the initial condition but vanishes as t — cxd. You can check 
for yourself by direct substitution that Eq. (50) is the solution to the inhomogeneous 
Heat Problem (45), i.e. the PDE, BCs and IC. 
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7.3 Heat sources 



7.3.1 Derivation 

Ref: Guenther & Lee p. 6 (Eq. 3-3) 

To add a heat source to the derivation of the Heat Equation, we modify the energy 
balance equation to read, 

change of 

^ ^ ^ heat in from heat out from ^ heat generated 

left boundary right boundary in segment 

segment in time At 

Let Q {x, t) be the heat generated per unit time per unit volume at position x in the 
rod. Then the last term in the energy balance equation is just QAAxAt. Applying 
Fourier's Law (1) gives 

/ du\ I Ou 

cpAAxu {x, t + At)- cpAAxu {x, t) = At A -Kq— - At A -Kq 



\ dx J ^ V dx / ^ , 

+QAAxAt 

The last term is new; the others we had for the rod without sources. Dividing by 
AAxAt and rearranging yields 

^ (a:, t + At) -u{x,t) ^ Kq f {^^)x+Ax ~ \ _^ 

At cp I Ax j cp 

Taking the limit At, Ax gives the Heat Equation with a heat source, 

du_ Q 

dt - ^dx' + cp ^^^^ 
Introducing non-dimensional variables x = x/l, i = nt/P gives 

dt dx'' Kcp ^ ' 

Defining the dimensionless source term q — PQ/ (Kcp) and dropping tildes gives the 
dimensionless Heat Problem with a source. 
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7.3.2 Solution method 

Ref: Guenther & Lcc p. 147 - 149, Myint-U & Dcbnath §6.7 (exercises) 

The simplest case is that of a constant source q = q{x) in the rod. The Heat 
Problem becomes, 

Ut = u^x + qix), <x <1, (54) 
u(0,t) = bi, u{l,t)^b2, t>0, 
u{x,0) = f {x) , < X <1. 

If q (x) > 0, heat is generated at x in the rod; if q (x) < 0, heat is absorbed. 

The solution method is the same as that for inhomogeneous BCs: find the equi- 
librium solution ue (x) that satisfies the PDE and the BCs, 

= u'^ + q{x) , < X <1, 
Ue (0) = hi, Ue (1) = 62- 

Then let 

V {x, t) — u {x, t) — Ue (x) . 
Substituting u {x, t) = v {x,t) + ue (x) into (54) gives a problem for v {x, t), 

Vt = fix, <x <1, 
v{0,t) = 0, v{l,t) = 0, t>0, 
v{x,0) = f{x)—UE{x), < a; < 1. 

Note that v (x, t) satisfies the homogeneous Heat Equation (PDE) and homogeneous 
BCs, i.e. the basic Heat Problem. Solve the Heat Problem for v (x, t) and then obtain 
u {x, t) — V {x, t) + Ue (x). 

Note that things get complicated if the source is time-dependent - we won't see 
that in this course. 

7.4 Periodic boundary conditions 

Ref: Guenther & Lee p. 189-190; for alternate method, see Guethcr & Lee p. 149 
and then p. 147 
[Sept 26, 2006] 

Above, we solved the heat problem with inhomogeneous, but time-independent, 
BCs by using the steady-state. We now show how to solve the heat problem with inho- 
mogeneous, but time-varying, BCs. We consider the heat problem with an oscillatory 
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(and periodic) BC, 



Ut = Uxx, < X < 1 

u{0,t) = AcosLot, u{l,t)^0, t>0, (55) 
u{x,0) = f{x), < X <1. 

The physical meaning of the BC ^(0,^) = Acosut is that we keep changing, in a 
periodic fashion, the temperature at the end ,t = of the rod. 

We don't expect the solution to be independent of time as t ^ oo, since we're 
changing the temperature periodically at one end. However, we do expect that after 
an initial transient time, the solution will become periodic with angular frequency uu, 
i.e. 

u {x, t) = V (x, t) + A (x) cos {ujt + (f) (x)) 

where v {x,t) ^ as t — >• cx) is the transient, A (x) cos {ut + (p (x)) is what we call 
the quasi-steady state, A (x) and </> (x) are the amplitude and phase of the quasi 
steady state. To solve the problem, the goal is to first find A (x) and (f) (x), and then 
v{x,t), if necessary. Often we might not care about the transient state, if we are 
more interested in the solution after " long times" . 

7.4.1 Complexify the problem 

We use the notation Re {2;}, Im{2;} to denote the real and imaginary parts of a 
complex number z. Note that 

Re {4 = hz + z*), cos^ = Re{e*^} 
Im{4 = )-.{z-z*), sin^ = Im{e^^} 

where asterisks denote the complex conjugate ((x + iy)* = x — iy). Thus, we can 
write our quasi-steady solution in terms of complex exponentials, 

A {x) cos {ujt + (j) {x)) = Re {A (x) e^^^^^e''^*} = Re {U {x) e'"*} 

where, for convenience, we have replaced A (x) e**^^^^ with the complex function U (x). 
We do this because complex exponentials are much easier to work with than cos (cut) 
and sin(c<;t). Note that U (x) has magnitude A{x) = \U {x)\ and phase 0(x) = 
arctan . The phase <p {x) delays the effects of what is happening at the end 

of the rod: if the end is heated at time t = ti, the effect is not felt at the center until 
a later time t — (j) (1/2) /u; + ti. The following result will be useful. 
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Lemma [Zero sum of complex exponentials] If, for two complex constants a, b, 
we have 

ae'""' + be-'"'' = (56) 

for all times t in some open interval, then a — b — 0. 
Proof: Differentiate (56) in time t, 

iuj {ae^^ - be-'""^) = 0. (57) 

Adding (56) to l/ia;x(57) gives 

2ae^'^* = 0. 

Since e*'^* is never zero (|e"^*| = 1), then a — 0. Prom (56), &e~"^* = and hence 
6 = 0. 

Note that we could also use the Wronskian to show this: 



W [e'^\ e-'^'] = det 



-2icu 7^ 



and hence e"^* and e ''^^ are linearly independent, meaning that a — b — 0. 



7.4.2 ODE and ICs for quasi-steady state 

Step 1. Pind the quasi-steady state solution to the PDE and BCs of the Heat Problem 
(55) of the form 

uss {x, t)^Re{U {x) e'^*} ^\{U {x) e^'"* + U* {x) e"*'"*) = A (x) cos {cut + (p{x)) 

(58) 

where U (x) is a complex valued function. Substituting (58) for u {x, t) into the PDE 
in (55) gives 

^ {iuU (x) e^'^* - iuU* (x) e-^'"*) = ^ {U" (x) e''^' + U"* (x) e"'^*) , 

for < X < 1 and t > 0. Multiplying both sides by 2 and re-grouping terms yields 

{iujU{x)-U"{x))e''^^ + {-iujU*{x)-U"*{x))e-''^' = 0, < a; < 1, t > 0. 

(59) 

Applying the Lemma to (59) gives 

iuU (x) - U" (x) = = -iujU* {x) - U"* {x) (60) 
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Note that the left and right hand sides are the complex conjugates of one another, 
and hence they both say the same thing (so from now on we'll write one or the other). 
Substituting (58) into the BCs in (55) gives 

I {U (0) e^"^* + U* (0) e"^"^*) = ^ (e^'"* + e"^'"*) , J {U (1) e^'"* + U* (1) e"^'"*) = 0, 

(61) 

for i > 0. Grouping the coefficients of e='=*'^* and applying the Lemma yields 

C/(0)=A [/(1) = 0. (62) 

To summarize, the problem for the complex amplitude U [x) of the quasi-steady- 
state uss {x,t) is, from (60) and (62), 

U" {x) - iuU (x) = 0; U (0) = A, U (1) = 0. (63) 

Note that (1 -I- i)^ = 2i, and hence 

iu;= ^{l + ifuj= + ^ 

Therefore, (63) can be rewritten as 

C/"-(y|(l + 0) U^O; C/(0)=A C/(1) = 0. (64) 

7.4.3 Solving for quasi-steady state 

Solving the ODE (64) gives 

U = ci exp (^-^^ (1 + + C2 exp (1 + i) x^ (65) 

where Ci, C2 are integration constants. Imposing the BCs gives 
^ = C/(0) =ci + C2, = [/ (1) = ci exp (-^f (1 + ^)) +C2exp(^y|(l + ?)^ . 

Solving this set of linear equations for the unknowns (ci,C2) and substituting these 
back into gives 

Aexp(yf (1 + i)) 



C2 = A - Cl 



exp(^(l + i)) -exp(-^(l + i)) 

Aexp (-/f (1 + 0) 



exp(yf (1 + i)) -exp(-^(l + i)) 
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Figure 4: At left, the magnitude of U{x) (solid) and (dash). At right, the 

phase of U (x). 

Substituting these into (65) gives 

ir _ ^4 ^^P (a/I + ^) - ^)) - (-/f + ^) (1 - ^)) 

exp (1 + i)) - exp (-^ (1 + i)) 
Therefore, the quasi-steady-state solution to the heat problem is 



( exp (y/f (1 + i) (1 - X)) - exp (- yf(l + .) (1 - X)) 

It is easy to check that uss {x,t) satisfies the PDE and BCs in (55). Also, note that 
the square of the magnitude of the denominator is 

exp (1 + 0^ ~ (^~'^^^^ + 0^ =2 ^cosh \/2a; - cos \/2uj^ > 



which is greater than zero since a; > and hence cosh v2ci; > 1 > cos v2ci;. In Figure 
4, the magnitude and phase of U (x) are plotted as solid lines. The straight dashed 
line is drawn with \U {x) \ for comparison, illustrating that \U {x) \ is nearly linear in x. 
The phase of U (x) is negative, indicating a delay between what happens at a point 
X on the rod and what happens at the end x — 0. 

7.4.4 Solving for the transient 

Step 2. Solve for the transient, defined as before, 

V {x, t) — u (x, t) — Uss {x-i t) . (66) 
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Substituting (66) into the heat problem (55), given that uss{x,t) satisfies the PDE 
and BCs in (55), gives the following problem for v (x, t), 



exp (1 + i)) - exp (- (1 + t)) 



Vt = Vxx-, < X < 1 

v{Q,t) = 0, 'u(l,i) = 0, i > 0, (67) 
v{x,{)) = /2(x), 0<x<l, 

where the initial condition /2 {x) is given by 
/a {x) = u (x, 0) - uss {x, 0) 

f{x)- Re |^^(^ + ~ ~ ^""P (- Vf (1 + (1 - ^)) 

The problem for v (x, t) is the familiar basic Heat Problem whose solution is given by 
(25), (24) with / (x) replaced by /s (x), 

oo „1 

V (x, t) = i?n sin {mrx) e~" *, Bn — 2 f2 {x) sin (nvra;) dx. 
n=l -^0 

7.4.5 Full solution 

The full solution to the problem is 

. exp ( vf (1 + 0(1 - X)) - exp (-71 (1 + .) (1 - .)) ^^^^ 

^ ' ^ exp(yf (l + i))-exp(-vf (1 + i)) 



+ ^5„sin {mrx) e""'"'* 

n=l 



The first term is the quasi-steady state, whose amplitude at each x is constant, plus 
a transient part v {x,t) that decays exponentially as t — >• oo. If the IC f {x) were 
given, then we could compute the B^s. 



7.4.6 Similar problem: heating/cooling of earth's surface 

Consider a vertical column in the earth's crust that is cooled in the winter and heated 
in the summer, at the surface. We take the x-coordinate to be pointing vertically 
downward with x = corresponding to the earth's surface. For simplicity, we model 
the column of earth by the semi-infinite line < a; < oo. We crudely model the 
heating and cooling at the surface as u{0,t) = Acosujt, where u = 2tx/t and the 
(scaled) period r corresponds to 1 year. Under our scaling, r = k x (1 year)//^. The 
boundary condition as a; — > oo is that the temperature u is bounded ("oo" is at the 
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bottom of the earth's crust, still far away from the core, whose effects are neglected). 
What is the quasi-steady state? 

The quasi-steady state satisfies the Heat Equation and the BCs, 

{uss)t = (uss)^^ , <x <oo (68) 
uss (0, t) — To + Ti coscut, uss bounded as x — > oo, t > 0. 

Wc use superposition: u {x, t) = Uq + Ui, where 

{uo)^ = (uo)^^ , {ui\ = (ui)^^ , < X < oo, i > 0. 

Uq (0, t) — To, Ui (0, t) — Ti coscut, Uq, Ui bouudcd as X — > oo. 

Obviously, Uq {x, t) = Tq works, and by uniqueness, we know this is the only solution 
for uo {x, t). To solve for ui, we proceed as before and let ui {x, t) = Re {U {x) e*'^*} 
to obtain 

U" (x) -iujU{x) ^ 0, < X < oo (69) 
U (0) = Ti, U bonded as x ^ oo, t > 0. 

The general solution to the ODE (69) is 

U = ci exp (^-^^ (1 + + C2 exp (1 + i) x^ . 

The boundedness criterion gives C2 = 0, since that term blows up as a; — > oo. The 
BC at the surface {x — 0) gives ci = Ti. Hence 

C/^Tiexp (^-^{l + t)x 

Putting things together, we have 

uss(x,t) = To + Re!^Tiexp(^-^^{l + i)x]e'' 



To + Tie'Vt^Re jexp (^-i'^x + iut 



= To + Tie"^"^ cos Y 1^ + j • (70) 

Uss {x,t) is plotted at various dimensionless times ut/Tr = 0, 1/4, 1/2,3/4, 1 in Figure 
5. Dashed lines give the amplitude Tq ± Tie~V^^ of the quasi- steady-state ugs {x,t). 

Physical questions: What is the ideal depth for a wine cellar? We want the wine 
to be relatively cool compared to the summer temperature and relatively warm to the 
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winter temperature, and yet we want the cellar close to the surface (to avoid climbing 
too many stairs) . Then we must find the smallest depth x such that the temperature 
ussixjt) will be opposite in phase to the surface temperature ^55(0,^). We take 
K, — 2x 10~^ cm^/s and I — 1 m. Recall that the period is 1 year, r = («/^^) (1 year) 
and 1 year is 3.15 x 10^ s. Prom the solution (70), the phase of uss {x,t) is reversed 
when 




Solving for x gives 

X = 7l\/2/u! 

Returning to dimensional coordinates, we have 

x' = lx = Itt\ = ^/^ui(YyeaI) = -^tt x (2 x IQ-^ cm^) x 3.15 x 10^ = 4.45 m. 
V 27r 

At this depth, the amplitude of temperature variation is 

Tie"^"^ = Tie"'' ^ OMTi 

Thus, the temperature variations are only 4% of what they are at the surface. And 
being out of phase with the surface temperature, the temperature at x' — 4.45 m is 
cold in the summer and warm in the winter. This is the ideal depth of a wine cellar. 

Note: for a different solution to this problem using Laplace transforms (not covered 
in this course), see Myint-U & Debnath Example 11.10.5. 

8 Linearity, Homogeneity, and Superposition 

Ref: Myint-U & Debnath §1.1, 1.3, 1.4 
[Sept 28, 2006] 

Definition Linear space: A set V is a hnear space if, for any two elements^ Vi, 
V2 &V and any scalar (i.e. number) A; e R, the terms V1+V2 and kvi are also elements 
of V. 

E.g. Let S denote the set of functions that are (twice-continuously differen- 
tiable) in x and (continuously differentiable) in t, for x E [0, 1], i > 0. We write <S 
as 

5 = {/ {x, t) I fxx, ft continuous for x e [0,1] , t >0} (71) 

If /i) /2 £ <S, i.e. /i, /2 are functions that have continuous second derivatives in space 
and continuous derivatives in time, and A; e R, then /i + /2 and kfi also have the 

"^Note that the symbol G means "an element of. So x G [0,1] means x is in the interval [0,1]- 
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same, and hence are both elements of <S. Thus, by definition, S forms a hnear space 
over the real numbers M. For instance, /i {x, t) — sin (ttx) e"'^^* and f2 — x'^ + t^. 
Definition Linear operator: An operator £ : V ^ W is linear if 

Civi + V2) = C{vi)+C{v2), (72) 
C{kvi) = kC{vi), (73) 

for all Vi, ^2 e V, /c e M. The first property is the summation property; the second is 
the scalar multiplication property. 

The term "operator" is very general. It could be a linear transformation of vectors 
in a vector space, or the derivative operation d/dx acting on functions. Differential 
operators are just operators that contain derivatives. So d/dx and d/dx + d/dy are 
differential operators. These operators are functions that act on other functions. For 
example, the x-derivative operation is linear on a space of functions, even though the 
functions the derivative acts on might not be linear in x: 

^{cosx + x^)^^{cosx) + l{x^). 

E.g. the identity operator I, which maps each element to itself, i.e. for all elements 
f G V we have I (f ) = v. Check for yourself that I satisfies (72) and (73) for an 
arbitrary linear space V. 

E.g. consider the partial derivative operator acting on S (defined in (71)), 

^ , du 

PM = ^, ueS. 

From the well-known properties of the derivative, namely 

d , . du dv ^ 
(^u + v) = — + u,veS, 
dx dx dx 

d dvj 

— (ku) = k—, ueS,keR 
dx dx 

it follows that the operator T> satisfies properties (72) and (73) and is therefore a 
linear operator. 

E.g. Consider the operator that defines the Heat Equation 

^ d 92 



dt 

and acts over the linear space S of functions C"^ in x, in t. The Heat Equation can 
be written as 
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For any two functions u,v E S and a real A; e R, 

Thus, L satisfies properties (72) and (73) and is therefore a hnear operator. 
8.1 Linear and homogeneous PDE, BC, IC 

Consider a differential operator jC and operators Bi, B2, T that define the following 
problem: 



PDE 
BC 
IC 



C{u) = h (x, t) , with 

Si(M(0,t)) = (7i(t), B2{u{l,t))^g2{t) 
l{u{x,0))^f{x). 



Definition The PDE is linear if the operator C is linear. 
Definition The PDE is homogeneous if h{x,t) = 0. 

Definition The BCs are hnear if Bi, B2 are linear. The BCs are homogeneous if 

91 (t) = 92 (t) = 0. 

Definition The IC is linear if T is a linear. The IC is homogeneous if / (x) ~ 0. 
8.2 The Principle of Superposition 

The Principle of Superposition has three parts, all which follow from the definition 

of a linear operator: 

(1) If £{u) — is a linear, homogeneous PDE and ui, U2 are solutions, then 
Ci'Ui + C2U2 is also a solution, for all Ci, C2 G M. 

(2) If ui is a solution to C{ui) = fi, U2 is a solution to C{u2) = f2, and C is 
linear, then ciUi + C2U2 is a solution to 

C (u) = ci/i + C2/2, 

for all ci, C2 e R. 

(3) If a PDE, its BCs, and its IC are all linear, then solutions can be superposed 
in a manner similar to (1) and (2). 



37 



The proofs of (l)-(3) follow from the definitions above. For example, the proof of 
(1) is 



£ (ciMi + C2M2) = £ (ciMi) + £ (C2M2) 

= c\L + C2L {U2) = ci • + C2 • = 

The first step follows from the addition rule of hnear operators; the second from the 
scalar multiplication rule, and the third from the fact that Ui and U2 are solutions of 
jC (u) — 0. The proofs of parts (2) and (3) are similar. 
E.g. Suppose ui is a solution to 

C (m) = Ut- u^x = 

Bi{u{0,t)) = u{0,t)=0 

B2{u{l,t)) = u{l,t)^0 

l{u{x,0)) = u{x, 0)^100 

and U2 is a solution to 

Bi{u{0,t)) = u{0,t)^100 

B2iu{l,t)) = u{l,t)^0 

X{u{x,0)) = u{x,0) = 

Then 2ui — U2 would solve 

Bi{u{0,t)) = u{0,t) = -100 
B2iu{l,t)) = u{l,t) = 
I {u {x, 0)) = u {x, 0) = 200 

To solve jC{u) — g with inhomogeneous BCs using superposition, we can solve 
two simpler problems: jC{u) — g with homogeneous BCs and jC{u) — with inhomo- 
geneous BCs. 

8.3 Application to the solution of the Heat Equation 

Recall that we showed that Un {x, t) satisfied the PDE and BCs. Since the PDE (8) 
and BCs (10) are linear and homogeneous, then we apply the Principle of Superposi- 
tion (1) repeatedly to find that the infinite sum u {x, t) given in (80) is also a solution, 
provided of course that it can be differentiated term-by-term. 
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9 Uniform convergence, differentiation and inte- 
gration of an infinite series 

Ref: Guenther & Lee Ch 3 talk about these issue for the Fourier Series, in particular 
uniform convergence p. 50, differentiation and integration of an infinite series on p. 
58 and 59, and the Weirstrauss M-Test on p. 59. 

Ref: Myint-U & Debnath §5.10, 5.13 (problem 12 p. 136) 

The infinite series solution (25) to the Heat Equation (8) only makes sense if it 
converges uniformly^ on the interval [0, 1] . The reason is that to satisfy the PDE, 
we must be able to integrate and differentiate the infinite series term-by-term. This 
can only be done if the infinite series AND its derivatives converge uniformly. The 
following results^ dictate when we can differentiate and integrate an infinite series 
term-by-term. 

Definition [Uniform convergence of a series] The series 

oc 
n=l 

of functions /„ (x) defined on some interval [a, b] converges if for every £ > 0, there 
exits an Nq {e) > 1 such that 



oo 

^ ^ fn {x , 

n=N 



< e, for all N > Nq. 
Theorem [Term-by-term differentiation] If, on an interval x e [(i,b], 

oo 

1. f (x) = fn {x) converges uniformly. 



n=l 



2. (x) converges uniformly, 

n=l 

3. {x) are continuous, 
then the series may be differentiated term-by-term, 

oo 



n=l 



^The precise definitions are outlined in §9.4 (optional reading). 

^Proofs of these theorems [optional reading] can be found in any good text on real analysis [e.g. 
Rudin] - start by reading §9.4. 
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Theorem [Term-by-term integration] If, on an interval x e [a, h] , 

oo 

1. f (x) — /„ (x) converges uniformly, 

n=l 

2. /„ (x) are integrable, 
then the series may be integrated term-by-term, 

rb rb 



f{x)dx = ^ / /„ {x) dx. 

n=l 



Thus, uniform convergence of an infinite scries of functions is an important prop- 
erty. To check that an infinite series has this property, we use the following tests in 
succession (examples to follow): 

Theorem [Weirstrass M-Test for uniform convergence of a series of functions] 
Suppose {fn {x)}'^^^ is a sequence of functions defined on an interval [a, 6], and sup- 
pose 

\fnix)\<Mr, ixe[a,b], n = l,2,3,...). 

Then the series of functions Xl^i fn (^) converges uniformly on [a, b] if the series of 
numbers converges absolutely. 

Theorem [Ratio Test for convergence of a series of numbers] The series 'Y2'^=i 
converges absolutely if the ratio of successive terms is less than a constant r < 1, i.e. 



< r < 1, (74) 



for all n > > 1. 

Note: N is present to allow the first N — 1 terms in the series not to obey the 
ratio rule (74). 

Theorem [Convergence of an alternating series] Suppose 

1. |ao| > l^il > \a2\ > • • • 

2. a2„-i>0, a2n<0 (n = 1,2,3...) 

3. lim„^oo On = 

then the sum 

oo 
n=l 

converges. 

Note: This is not an absolute form of convergence, since the series Yl'^=i ("-'-)" /"^ 
converges, but Yl'^=i ^/^ does not. 
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9.1 Examples for the Ratio Test 

E.g. Consider the infinite series of numbers 



id 



n=l 



(75) 



Writing tliis as a series X^^i On, we identify = (1/2)"'. We form the ratio of 
successive elements in the series, 

|a„+i| (1/2)"+^ 1 ^ ^ 



Kl (1/2)" 2 

Thus, the infinite series (75) satisfies the requirements of the Ratio Test with r — 1/2, 
and hence (75) converges absolutely. 
E.g. Consider the infinite series 



OO -. 

E- 

^ n 

n=l 



(76) 



Writing this as a series X^^i '^^ identify an — l/n. We form the ratio of successive 
elements in the series, 

_ 1/ (n + 1) _ n 
|a„| 1/n n + 1 

Note that lim„^oo = Ij and hence there is no upper bound r < 1 that is greater 
than |a„+i| / |a„| for ALL n. So the Ratio Test fails, i.e. it gives no information. It 
turns out that this series diverges, i.e. the sum is infinite. 
E.g. Consider the infinite series 



n=l 



Again, writing this as a series J2'^=i ^n-i we identify a„ 1/n^. We form the ratio of 
successive elements in the series. 



|a„+i| l/(n+l)^ n 



2 



\an\ l/ri^ {n+lf 

2 

Note again that hm^^oo , ",^2 ~ 1, and hence there is no r < 1 that is greater than 
|a„+i| / \an\ for ALL n. So the Ratio Test gives no information again. However, it 
turns out that this series converges: 



n=l 
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So the fact that the Ratio Test fails does not imply anything about the convergence 
of the series! 

Note that the infinite series 

E - (78) 

n=l 

converges for p > 1 and diverges (is infinite) for p < 1. 

9.2 Examples of series of functions 

Note, 

oo 

sin (riTTx) 

rt=l 

does not converge at certain points in [0, 1], and hence it cannot converge uniformly 
on the interval. In particular, at a; = 1/2, we have 

oo oo 

sin (nTTx) ^ sin = l + 0- l + + l + 0- H 

n=l n=l 

The partial sums (i.e. sums of the first n terms) change from 1 to forever. Thus the 
sum does not converge, since otherwise the more terms we add, the closer the sum 
must get to a single number. 
Consider 

n=l 

Since the n'th term is bounded in absolute value by and since Yl^=i 1/^^ 'con- 

verges (absolutely), then the Weirstrass M-Test says the sum (79) converges uniformly 
on [0, 1]. 

9.3 Application to the solution of the Heat Equation 

Rcf: Guenther & Lee p. 145 - 147 

We now use the Weirstrass M-Test and the Ratio Test to show that the infinite 
series solution (25) to the Heat Equation converges uniformly, 

oo oo 

u (x, t) = ^Un (x, t) = ^Bn sin (nTTx) e"""''"*, (80) 

n=l n=l 

for all positive time, i.e. t > to > 0, and space x e [0, 1], provided the initial condition 
/ (x) is piecewise continuous. 
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To apply the M-Test, we need bounds on |it„ (x, t)\, 



\un {x, t) \ — Bn sin {nnx) e 



for all X e [0, 1] . (81) 



We now need a bound on the Fourier coefficients Note that from Eq. (24), 



2 / sin (rmrx) f (x) dx 



<2 \sm {mnx) f {x)\dx < 2 \f {x)\dx, (82) 



\B„ 



for all X e [0, 1]. To obtain the inequality (82), we used the fact that |sin (mnx)] < 1 
and, for any integrable function h{x), 



rb rb 

/ h{x)dx < \h{x)\dx. 

J a J a 



You've seen this integral inequality, I hope, in past Calculus classes. We combine (81) 
and (82), to obtain 

\Un{x,t)\<Mn (83) 



where 



M„= (2 j^' \fix)\dx^ e 



To apply the Weirstrass M-Test, we first need to show that the infinite scries of 
numbers Xl^i ^« converges absolutely (we will use the Ratio Test). Forming the 
ratio of successive terms yields 



^(.„2_(„+,)2)^2,^ ^ g-(2n+l),r% < ^-.H, n = 1, 2, 3, ... 



Thus, by the Ratio Test with r = e " < 1, the sum 



n=l 



converges absolutely, and hence by Eq. (81) and the Weirstrass M-Test, (x, t) 

converges uniformly for a; G [0, 1] and t > to > 0. 

A similar argument holds for the convergence of the derivatives Ut and Uxx- Thus, 
for all i > to > 0, the infinite series (80) for u may be differentiated term-by-term 
and since each Un (x, t) satisfies the PDE and BCs, then so does u (x, t). Later, after 
considering properties of Fourier Series, we will show that u converges even at t = 
(given conditions on the initial condition f (x)). 
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9.4 Background [optional] 

[Note: You arc not responsible for the material in this subsection 9.4 - it is only added 
for completeness] 

Ref: Chapters 3 & 7 of "Principles of Mathematical Analysis" , W. Rudin, McGraw- 
Hill, 1976. 

Definition Convergence of a series of numbers: the series of real numbers 
converges if for every £ > 0, there is an integer N such that for any n> N, 



oo 
m=n 



E 



< s. 



Definition Absolute convergence of a series of numbers: the series of real numbers 
'^n is said to converge absolutely if the series |art| converges. 

Definition Uniform convergence of a sequence of functions: A sequence of func- 
tions {/„ (a;)}^-^ defined on a subset E C R converges uniformly on E to a function 
/ {x) if for every e > there is an integer N such that for any n > N, 

\fn{x)-f{x)\ < e 

for all X eE. 

Note: pointwise convergence does not imply uniform convergence, i.e. in the 
definition, for each e, one works for all x in E. For example, consider the sequence 
of functions {x"'}^^^ on the interval [0, 1] . These converge pointwise to the function 

0, 0<x<l, 

1, X = 1 

on [0, 1] , but do not converge uniformly. 

Definition Uniform convergence of a series of functions: A series of functions 
Yl'^=i fn (x) defined on a subset E C R converges uniformly on E to a function g (x) 
if the partial sums 

m 

Sm {x) = ^ /n (x) 
n=l 

converge uniformly to g (x) on E. 
Consider the following sum, 

oo . X 

Ecos [nirx) 
n 

n=l 

docs converge, but this is a weaker form of convergence (pointwise, not uniform or 
absolute), by the Alternating Series Test (terms alternate in sign, absolute value of 
the terms goes to zero as n — > oo). Consult Rudin for the Alternating Series Test, if 
desired. 
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